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Abstract 

One of the simplest time-dependent solutions of M theory consists of nine-dimensional Euclidean 
space times 1+1-dimensional compactified Milne space-time. With a further modding out by Z2, 
the space-time represents two orbifoldplanes which collide and re-emerge, a process proposed as 
an explanation of the hot big bang Q,y,Q. When the two planes are near, the light states of the 
theory consist of winding M2-branes, describing fundamental strings in a particular ten-dimensional 
background. They suffer no blue-shift as the M theory dimension collapses, and their equations of 
motion are regular across the transition from big crunch to big bang. In this paper, we study the 
classical evolution of fundamental strings across the singularity in some detail. We also develop 
a simple semi-classical approximation to the quantum evolution which allows one to compute the 
quantum production of excitations on the string and implement it in a simplified example. 
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I. INTRODUCTION 



Recently, an M theory model of a big crunch/big bang transition has been proposed. 
The model consists of two empty, parallel orbifold planes colliding and re-emerging. Close 
to the collision, the spectrum of the theory splits into two types of states. The light states 
are winding M2-branes describing perturbative string theory, including gravity. The massive 
states, corresponding to higher Kaluza-Klein modes associated with the M theory dimension, 
are non-perturbative DO-brane states whose masses diverge as the collision approaches. In 
Ref . j^f , it was shown that the equations of motion of the winding M2-branes are regular at 
the brane collision and yield an an unambiguous classical evolution across it. Since these 
states describe gravity in quantized string theory, the suggestion is that gravity is better- 
behaved in this situation than it is in general relativity. If this is true then, rather than 
seeding a gravitational instability, the massive Kaluza-Klein states may simply decouple 
around the collision. Provided their density is negligible in the incoming state, as it is in the 
ekpyrotic and cyclic models , then the entire transition may be describeable using 

the perturbative string theory states alone. 

In this paper, we will assume the massive states can be neglected and focus on the 
evolution of the fundamental string states across a big crunch/big bang transition of this 
type. By solving the string evolution equations numerically we show how the higher string 
modes become excited in a well-defined way as strings cross the transition. We also develop 
a simple approximation whereby quantum production of excited modes on the string may 
be computed in a semi-classical manner. Related work may be found in Ref. Earlier 
work on strings in cosmological backgrounds is reviewed in Ref. 0. 

The background space-time we are interested in is most simply understood in eleven- 
dimensional terms. Near the orbifold plane collision, the line element reduces to that for a 
Z 2 -compactified Milne universe times R 9 , namely 

9 

ds 2 = -dt 2 + t 2 d9 2 + E( rf ^) 2 , (1) 

where < 9 < 9 , and a Z 2 reflection is imposed about each end point. For t < (t > 0), the 
two orbifold planes approach (recede) with relative rapidity 9q. Away from the singularity 
at t — the space-time is flat and hence an automatic solution of any theory governed by 
field equations involving purely geometrical terms 0. At t — the metric degenerates 
and the equations of general relativity become singular. Nevertheless, one can analytically 
continue the background solution through t = 0. As mentioned, the equations governing the 
winding M2-branes are regular at t — 0, so the M2-branes evolve smoothly across the big 
crunch/big bang transition. 
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The winding M2-branes reduce to fundamental strings when the M theory dimension 
is small, so one can also describe the situation in ten-dimensional terms. In string frame, 
the dimensionally-reduced metric takes the form g^ u = | t/t s 1 77^^ , with a time- dependent 
dilaton, <fi = |ln|£/i s |. Here is the string coupling and t s is by definition the time 
when the string coupling is unity [df. For \t\ < f S) stringy interactions are weak and one 
can, to a first approximation, treat the strings as free. For \t\ greater than or of order 
t s , however, the string theory is at strong coupling and one must switch to an eleven- 



dimensional supergravity description 10|,lll|- Recently, techniques for solving the relevant 



higher- dimensional Einstein equations have been developed 12j , which can be applied to this 
situation. As long as one is considering long-wavelength cosmological perturbations, evolving 
classically, the transition from the supergravity regime to the string theory (a' expansion) 
regime appears unproblematic: the eleven- dimensional Einstein equations should reduce to 
the appropriate ten dimensional Einstein-dilaton effective theory as the branes become near. 
In this paper, our main focus is on the classical dynamics of the string near t = 0. In a 
companion paper we provide a complementary treatment, studying the stringy (3 function 
equations in the usual a' expansion to see whether higher order corrections are significant 
when the string scale crosses the Hubble radius 3|. 

In the eleven dimensional picture, there is no dilaton field. The unique mass scale is set 
by the M2-brane tension /i 2 = Mf l} where M n is the eleven- dimensional Planck mass [14]. 
Winding membranes, of length 9o\t\ behave like strings with a time-dependent tension, /ii = 
Mii^ol^l, m Minkowski space-time. One can equally well view the strings as having a fixed 
tension, \i\ = Mf^oltsl, but living in a string-frame background metric, g s = \t/t s \rj^ u . The 
string coupling constant is set by the size of the M theory dimension in eleven- dimensional 
Planck units: = (M n 9o\t\)% = \t/t s \l. 

Starting well away from the collision, at large times \t\ the system is described by eleven- 
dimensional supergravity. As one approaches the collision the string coupling falls below 
unity, when \t\ falls below 

t s = e^ l M u l . (2) 

When \t\ = t s , the string mass scale is M n and the string length is l s ~ ~ 6* t s . For 
small #0, l s is far smaller than the Hubble horizon scale t s so the characteristic oscillations 
of the string are little affected by the expansion of the universe. Corrections due to the 
background space-time curvature are small and the usual a' expansion holds good. However, 
as \t\ decreases further, the physical size of an oscillating string remains approximately 
constant while the physical Hubble radius falls as |t|(|t|/t s )^. The string length crosses the 
Hubble radius at a time 

t x ~elt s = e~ li M u \ (3) 
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FIG. 1: Three phases of an orbifold plane collision, with rapidity 0q < 1. Cosmological perturba- 
tions are created when the effective theory is the low energy limit of M-theory in eleven dimensions. 
After the orbifolds get closer and the string coupling (g s ~ |i/i s | 3 / 2 ) falls below unity, the system 
is better described in terms of perturbative string theory and the usual expansion in a'. A second 
transition occurs when the string scale leaves the Hubble radius. Here the string tension becomes 
dynamically subdominant and the theory is better described by an expansion in 1/a! . In this 
regime the string is better-described as a collection of independent 'bits' whose evolution across 
the singularity is smooth. After t = the string tension is reconstituted and the system runs back 
through the two abovementioned transitions. 



Thereafter, the string tends towards an ultralocal evolution in which the curvature of the 

string is unimportant and each bit of the string evolves independently. This is the regime of 

the 1/a' expansion discussed in Ref. 3]. Figure 1 illustrates the three regimes: from strong to 

_i 

weak coupling at t s ~ ^ 1 M{] 1 and from flat space to ultralocal evolution at tx ~ @o 3 . 
For small 9q, as we shall assume throughout this paper, these times are well-separated. 

In this paper, we wish to study the weak coupling regime, \t\ < t s , during which the 
strings cross the singularity. In this regime, for \t\ < tx, the usual a' expansion fails and we 
must replace it with an expansion in the string tension, proportional to 1/a'. Our ultimate 
goal is to compute a 'mini S-matrix', evolving weakly-coupled incoming string states from 
an initial time just after — t s to a final time just before +t s . In this paper we shall make 
modest progress towards this goal by computing the 'mini S-matrix' in a first, classical 
approximation. 

Any classical solution of the string equations also solves the Heisenberg operator equations 
for the string field x M (r, a), to leading order in h. For highly excited string states, involving 
large occupation numbers, the classical approximation should be reasonable. However, to 
describe gravity we are really more interested in the low-lying states. For these states, 
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quantum corrections are likely to be significant and there is less we can say with precision. 

Nevertheless, there are arguments which suggest that quantum corrections may be man- 
ageable in the relevant regime. The usual a' expansion is based on the approximating 
space-time as locally flat, and solving the string equations of motion in an expansion in 
the space-time curvature. This approximation clearly fails to describe even the classical 
motion of the fundamental strings near t = 0. Hence the failure of the a' expansion in the 
quantum theory does not, on its own, imply that quantum corrections to the string evolu- 
tion are large [7]. Second, the key phenomena relevant to the small time regime take place 
when the background curvature length (i.e. the Hubble radius) falls below the string scale. 
The string dynamics becomes ultralocal in this regime: in cosmological terms, the string is 
'super-horizon' and hence may be expected to evolve in a classical manner. 

In flat spacetime, the most dramatic consequence of quantizing the string is that the first 
excited states of the string are massless instead of massive as a classical treatment would 
suggest. This is due to quantum mechanical renormalization of the mass squared operator, 
m 2 = (p ) 2 — p 2 , where p M is the string's centre of mass momentum. The renormalized mass 
squared operator is lowered by one quantized unit relative to the normal ordered classical 
expression, and this lowering causes the first excited states to become massless. As a result, 
the string centre of mass trajectory x^ m = + p^r is altered from the classical expression 
by (p ) 2 — > (p ) 2 — M 2 with M a constant so that the quantum corrected trajectory is 
null instead of spacelike. Note, however, that this shift in the energy p° merely alters 
the time coordinate of the centre of mass of the string and it does not affect the string 
spatial coordinates and momenta, which obey exactly the same equations as in the classical 
theory. Again, this discussion suggests that quantum effects may be relatively modest in 
the situation we are interested in. 

In this paper, we focus on the classical evolution of the string across t = 0. Starting at 
a time t of order —t 8 , we expect the string to oscillate nearly adiabatically and with almost 
fixed physical size until the string scale crosses the Hubble radius at time t = —tx- For 
small #0; t s 3> tx, hence there is a large range of time over which adiabatic evolution holds. 
This allows us to clearly identify the incoming and outgoing states in terms of the usual flat 
space-time modes. As we have already mentioned, the natural size scale for the incoming 
states we are interested in (i.e. the graviton, dilaton or antisymmetric tensor states) is the 
string scale at t s , 

Ax~ M u l ~6 t s . (4) 

By the same token, the natural measure for the spatial momentum of the string state at this 
time is 

P ~ M u ~ (0 O O -1 > (5) 
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and we shall always measure the momentum of the incoming string in this basic units. 

One of the nice features of the problem at hand is that the string tends to simple flat 
space evolution at large asymptotic times: 

^( rj<T )-> ( in V(r,a) *->-*., (6) 

and 

XA»( T)<7 )_ (°«*V(T )( 7) (7) 

where m x fJ- (r, a) and ott *x At (r, a) are free string fields evolving adiabatically in the usual flat 
space string modes. The limits should, strictly speaking, be expressed in terms of a t c < t s 
so that we are always studying times for which the string coupling is smaller than unity. 
In practice this will not be an important distinction since as long as 9q is small, the string 
states we study tend to their asymptotic behavior at times t x <C t s - 

In a linear field theory in a time-dependent background, where the field tends to the 
usual Minkowski evolution in the asymptotic past, all of the information regarding quantum 
amp litudes may be obtained from real solutions of the classical field equations (see e.g. 
[la]). One simply expresses the incoming field as a linear combination of creation and 
annihilation operators, defined by their action on the incoming vacuum state, multiplied 
by the appropriate positive and negative frequency modes. The classical field equations 
then determine the quantum Heisenberg field for all time. Using the time-dependent field, 
one can then compute any correlation function of interest, evaluated in any incoming state 
expressed in terms of incoming creation operators acting on the incoming vacuum. 

The main difference in our situation is that the string evolution is nonlinear. Nevertheless, 
the nonlinearity turns out to be of a very simple form, where one can apply a semi-classical 
approximation in a self-consistent way. One can follow a similar procedure to that in the 
linear case to compute the quantum production of string excitations due to passage across 
the singularity. We shall explain this calculation, for a specially simple circular loop 

- in the final section of this paper. The calculation may be generalized to higher modes of 
the string and we shall do so in a future publication. 

The outline of this paper is as follows. In section II we present the equation governing 
winding membranes, discuss their adiabatic solutions and define various quantities of in- 
terest. In section III we study the classical evolution of a circular loop, representing the 
dilaton and its massive counterparts. In section IV we study the classical evolution of a 
rotor, representing the modes of maximal spin (for example, the graviton), at each mass 
level. Section V describes the classical evolution of a state of intermediate angular mo- 
mentum, representing an antisymmetric tensor state and massive analogs. Section VI gives 
an example of a quantum transmutation amplitude computed for an incoming dilaton-like 
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state. In section VII we give a brief summary with future proposals. In Appendix A we give 
an analytic solution method for the classical string equations, representing an expansion in 
powers of the string tension, i.e., inverse powers of a', which accurately describes the passage 
of a circular loop across the singularity. Intriguingly, in this regime polylogarithm functions 
enter, perhaps hinting at some deeper analytic solution yet to be found. We also present an 
analytic solution to second order in 1/a' for a more general string configuration. Finally, 
Appendix B is devoted to some numerical checks on the right/left mover decomposition for 
the special case of a classical rotor. 

II. EQUATIONS OF MOTION FOR A WINDING MEMBRANE 

Our starting point is the Nambu-Goto action for a membrane, 



where a a (a = 0, 1, 2) are the three world-volume coordinates, G a b = d a x A dbX E ' qab is the 
induced metric on the world- volume, x A = (t, x, y) are the space-time embedding coordinates 
and fi2 is the membrane tension. For a winding membrane in its lowest Kaluza-Klein state, 
we can set a 2 = 9 in the Milne metric The action (jSJ) then reduces to 



where now G^l is the induced metric on the two-dimensional world sheet with coordinates 
cr° = r and a 1 = a. Henceforth dots shall denote derivatives with respect to r and primes 
derivatives with respect to a . 

The action (JHJ) may be viewed as describing a string moving in flat space-time with a 
tension /i2#o|£| tending to zero as the string approaches the singularity. Hence one expects 
all points on the string to move with the speed of light in this limit. Equivalently, the 
same action may be viewed as describing a string of fixed tension ^^ots, moving in a time- 
dependent background with metric = (i/isl^. We shall adopt this latter point of view 
throughout the paper. 

The action @ is invariant under reparametrizations of the string world sheet. In partic- 
ular, we can choose timelike gauge x° = t — r, and also x ■ x ' = 0. In this gauge, the string 
spatial coordinate x(t, a) obeys the following classical equations of motion: 




(8) 




(9) 




(10) 
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where e is an auxiliary quantity (roughly speaking, the 'relativistic energy density' for the 
string), defined by 



W (11) 



X 2 ) 



Equations (fTUJl and (fTTjl imply that 



d t {e*) = \{l-x*)e\ (12) 

which will be useful later. Notice in particular that as t tends to zero, from (jlOj) e tends to 
a constant. Hence from from (JT2)l the speed of the string tends to unity. As pointed out in 
Ref. p|, for generic string states, equations (|1U|) are regular for all t. 

As in flat space-time, it is useful to rewrite these equations in terms of left and right 
moving modes, defined by 

f=x— — x', l = x + —x". (13) 

It is easy to show that I and r are unit vectors, r 2 = 1 = Z 2 , as in flat space 0]. For a 
closed string, they each describe a closed curve on a unit sphere. Notice also that whereas 
the timelike gauge we have chosen is invariant under reparameterizations of a, the the left 
and right movers are themselves reparameterization invariant. Many of the properties of 
oscillating loops can be seen most directly by picturing the left and right movers as curves 
on a unit sphere. In three spatial dimensions, even in flat space-time such curves generically 
cross. Where this happens the left and right movers coincide, and it follows that the string 



moves at the speed of light for an instant 18]. As we have already mentioned, in the 



background of interest, at t = every point on the string must move at the speed of light. 
Hence, at this moment, the the left and right mover curves must actually coincide. 

In terms of the left and right movers, the equations of motion (jl(J|) read 

?+ —r" + (j- (f ■ f)rj = 0, 

l-^l' + ^ t (r-(f.l).l) = 0. (14) 

In the limit t — > 0, the last terms in each equation force 

r^l, (15) 

so that, as explained earlier, all points on the string reach the speed of light. Conversely, 
in the limit of large times, when a loop is well inside the Hubble horizon, one expects it to 
evolve as in flat space-time. As explained below, when we convert to proper time and space 
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coordinates, the right and left movers defined in (JT3j) correspond to flat space-time right and 
left movers after a reparameterization of o. Our numerical calculations verify that f and 
I tend towards fixed curves at large times, providing a straightforward matching onto flat 
space-time solutions in this asymptotic regime. 

The equations of motion are non-linear and hard to solve analytically. We have therefore 
resorted to a numerical study of a variety of cases. We have also developed analytic approx- 
imations in the large time and small time limits, which we compare to the numerical results. 
Recall that the string length is ~ M^ 1 when we start our evolution, and the Hubble radius 
is larger by a factor 9q . As time runs forward, the comoving Hubble radius shrinks as \t\ 
and the comoving loop size grows as \t/t s \~2. Once the string crosses the Hubble radius, its 
vibrations cease and the string follows a kinematical super-Hubble evolution. 

The 'freezing' of the string outside the Hubble radius is, as we will discuss below, like 
a sort of measurement process. Depending on the phase of oscillation of the string, either 
the string coordinate can be frozen, or its momentum. In either case, the conjugate variable 
then acquires a large kick following t = 0. After crossing the singularity, the loop re-enters 
the Hubble radius and starts oscillating as in flat space once again (see Figure HJ). 



A. Asymptotic states: fiat space description 

At early and late times the string loop is well inside the Hubble radius and we expect 
it to follow standard flat space-time evolution, which we now pause to review. In standard 
flat coordinates, 

ds 2 fl = -dT 2 + (dX) 2 , (16) 

we can choose timelike, orthonormal gauge, T = r and (8tX) ■ (d a X) = 0. The string then 
evolves according to the wave equation 

d 2 X = d 2 J, (17) 

with the constraint 

(d T X) 2 + (d a X) 2 = 1. (18) 
The usual left and right movers are given by 

R = d T X - d a X, L = d T X + d a X, (19) 

and the equations of motion may be written 

d T R + d a R = 0, d T L - d a L = 0. (20) 
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The solutions may be expressed as a sum over Fourier modes, 

+00 +00 
R\T-a)= <e" in(T " ff) , L\T + a)= £ < e - m(T+<T) , (21) 

n=— 00 n=— 00 

where reality imposes (a l n )* = a l _ n and (a l n )* = a l _ n . To quantize the string, these Fourier 
parameters are promoted to operators. After following the canonical procedure, and reg- 
ularizing and renormalizing the nonlinear constraints, one finds there are three massless 
excitations, consisting of states of the form a^a^lO) where |0) is the oscillator vacuum 
state. These states consist of a space-time scalar (the dilaton), a symmetric and traceless 
tensor (the graviton) and an antisymmetric tensor. These modes are the most relevant cos- 
mological perturbation theory and evolving perturbations through the bounce. Hence, the 
rest of the paper will be dedicated to their classical analogs. We shall not, however, study 
the important issue of mass renormalization, which we defer to future work. 

The classical string configuration with no angular momentum, hence corresponding to a 
dilaton-like state, consists of a circular loop. In its rest frame, a loop in the XY plane takes 
the form: 

X D = cos(T)(cos((x),sin(a),0,...,0), < a < 2tt. (22) 

In contrast, the state with maximal angular momentum for a given energy, analogous to the 
graviton, takes the form of a rotor, a spinning doubled line: 

X G = cos((x)(cos(r),sin(T),0, . . .,0), < a < 2vr. (23) 

The left and right movers are easily calculated: they are circles in the XY plane which are 
parallel in the case of the dilaton-like state and antiparallel in the graviton-like state. A 
state of intermediate angular momentum may be constructed by taking the left and right 
movers to trace out two circles in perpendicular planes. Choosing the right mover in the 
XY plane, and the left mover in the XZ plane, the classical solution is 

X A = i(2cos(T)sin((x), cos(T-a), cos(T + a), 0, 0), (24) 

The three solutions and their description in terms of left and right movers on the sphere are 
depicted in Figure El 

In order to connect these flat space-time solutions to those solutions in an expanding 
universe, we need to relate comoving coordinates x and conformal time t to physical coor- 
dinates X and proper time T. Since dX = adx and dT = a dt, the factors of the scale 
factor a cancel out in the velocity. One can also show, using the result of Ref. ^| that the 
time-averaged velocity squared of a loop in its center of mass frame is |, that (|T2^) implies 
that J eda oc \t\2 oc a{t). Using this result, one can show using an appropriate reparam- 
eterization of a that the left and right movers (fTB*)) map precisely onto the flat space left 
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FIG. 2: Classical flat space-time solutions consisting of the lowest mode excitations with minimal, 
maximal and intermediate angular momentum. To get zero angular momentum the right and 
left movers lie in the same plane with parallel orientation. This configuration corresponds to 
an oscillating circular loop. The solution with maximal angular momentum is similar but has 
the left and right movers oppositely oriented. This solution describes a spinning doubled line, 
or rotor. A classical solution with intermediate angular momentum is achieved by taking the 
right mover circle to lie in a plane orthogonal to the left mover circle. The solution is three- 
dimensional. However, under a suitable projection it appears as an elliptical loop, spinning through 
an intermediate doubled-line configuration. The WV plane is spanned by V = (0, 1, 1,0, ...,0)/\/2 
and W = (0,-1,1,0, .. .,0)/y/2. 

and right movers This makes it easy to send in states of any desired asymptotic form, 
and to read off the states in which they come out. A further convenience of the timelike 
orthonormal gauge we use is that we do not need to actually perform the reparameterization 
of a. At any value of a one can read off the left and right movers as functions of the time t. 

B. Conventions and important quantities for the incoming/outcoming modes 

As noted above, we will generally choose a to run from to 2tt. Moreover, in our 
convention capital letters (like X, E, etc.) denote flat space-time quantities, i.e. those which 
correspond to flat space-time variables in the regime where the expansion of the universe is 
adiabatic (i.e. where the string loop is well inside the Hubble radius). On the other hand, 
lower case symbols (like x, e, t, etc) describe coordinates or quantities in the cosmological 
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background. To specify an incoming mode, for the circle or rotor we will assume the initial 
string configuration is in the xy plane, with some center of mass velocity in the z direction. 
We define the latter by 

Jezda 

v * = -JTdV' (25) 

calculated at a time where we start the evolution, which we shall formally denote —t s . The 
flat space solutions used to specify initial conditions are obtained by Lorentz boosting those 
given above. Moreover, we will define the comoving string size at t s to be 

Ax = — f da\x'\. (26) 

and we shall choose this to be unity in the initial state. For each choice of the initial time t s , 
there will be an associated orbifold rapidity 6q =t~ 1 . Finally, to get a quantitative measure 
of how much energy was produced (or lost) during the t — transition, one can define a 
flat-space-time energy 

E = J dae\t s t\- 1/2 (27) 

which approaches a constant for large values of \t\, as explained above. The energy ratio 
between an incoming mode and an outgoing one (rj = E in /E out ) provides a measure of how 
much energy was produced during the string's passage across the singularity. 

We will now consider separately the behavior of the circle, the rotor and the classical 
analog of the antisymmetric tensor. 



III. CIRCULAR LOOP 

The simplest nontrivial case to study is a circular loop, where the symmetry reduces 
the problem to a single dynamical variable. Employing the symmetric ansatz x(t, a) = 
p(r) ( cos(cr), sin(a), 0, ) where p is the comoving radius of the loop, in the gauge t = r 
the string action Q reduces to 

S = -2tt^ 2 / dt6 \t\ |p| \A ~ P 2 - ( 28 ) 

To minimize clutter, from now on we shall work in time and length units in which 2irp2do = 1- 
The canonical momentum p conjugate to p and the Hamiltonian are found respectively to 
be 

dp \t\S\p\ 
d 

and Hamilton's equations are 



dp _ _f_P dp_p_ . . 

dt~ H ' dt ~ H- m 
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FIG. 3: Numerical solution to equations (|29 |) - (j31|) . For large \t\, the comoving radius p oscillates, 
with an envelope growing as \t s /t\ l l 2 . The evolution across t = is smooth and, in this particular 
example, it results in a larger amplitude oscillation after the transition. 

Finally, because the Hamiltonian is expicitly time- dependent, we obtain 

dH 2 9 

-w = 2t " ■ (31) 



The qualitative properties of the solutions to (J2HJ) - d3H) are easily seen. The energy density 
e = H takes its minimal value at t = 0, where the radial speed dp/dt reaches unity. This last 
effect triggers interesting effects which we shall describe below. The solutions are regular 
for all t except for the special case where the momentum p vanishes at t — 0. For that case, 
(13 1|) implies that H ~ p \t\ at small t. From (}3T)j) and (|29| . one finds that p is very mildly 
non-analytic, p ~ — |t 2 sign(t), whereas p ~ p — \t 2 / po is regular. 

Equations (|29|) - (|31|) are invariant under the rescaling t — > At, e — * A 2 e, p — * A 2 p and 
p — > Ap, hence solutions for loops of different sizes are trivially related, and inequivalent 
solutions are labelled by only one parameter, which may, for example, be taken to be the 
asymptotic phase of the oscillation. 

The motion of the loop is simplest to describe in two asymptotic regimes, corresponding 
to large and small \t\, when the loop's radius p is well within, or outside, the Hubble 
radius. In the first regime the loop oscillates with fixed amplitude and period in proper 
time. In comoving coordinates and conformal time the oscillation amplitude changes as 
a(t) -1 = (ts/ltl) 1 / 2 , and the frequency changes as a{t) oc ( | ^ | /t s ) 1,/2 , as illustrated in Figure 
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FIG. 4: The radius (solid black line) and the velocity (blue dashed line) are plotted for two different 
circular loops. In the left plot the singularity is crossed after a maximum in the radius and almost at 
a minimum in the velocity. In this particular case, there is a loss of energy after the transition: the 
outgoing mode has a smaller amplitude than the incoming one. The right plot shows the opposite 
situation, with the radial coordinate approaching a maximum at the singularity, and energy being 
produced. 

The energy density e scales as a(t) = (Itl/t^) 1 / 2 in this regime. 

Once the co moving Hubble radius \t\ falls below the loop radius p, the loop enters a new 
kinematical regime in which the tension plays a subdominant role. As the time tends to 
zero, all points on the string approach the speed of light. The loop receives a 'kick' from 
passing through t = and emerges with a shift in its its oscillation phase and a net energy 
gain or loss. In the case of a circular loop, if the radius is expanding at t — the loop gains 
energy, whereas if it is contracting the loop loses energy. Time reversal, which is a symmetry 
of the equations, relates these two situations. Figure E| shows two examples. 

A. Asymptotic states at large \t\ 

In order to parameterize the incoming and outgoing states, it is helpful to perform a 
canonical transformation to new coordinates (0, V): 



V 




tan(0) 



p\t 



(32) 



2\t 



P 
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At large times, the angle represents the oscillation phase while the conjugate momentum 
V tends to a constant measuring the energy stored in the loop. The Hamiltonian equations 
now read 

dV V (OA,\ 

These equations are readily solved at large times, giving 
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V^V 00 = const., <f>^ J^±dt~±^J£- + ( j> , (34) 

where O is a constant phase. From it follows that the energy (|77jl E oc \/V and hence 
it tends to a constant at large times, as expected. In the usual comoving p coordinate, the 
asymptotic behavior of the solution at large times is 



W^ Sin (W^ + °J- (35) 

Upon identifying proper time and radius, this takes the form of the flat space solution (|22j). 
up to an arbitrary phase. Due to the rescaling symmetry discussed above, the only nontrivial 
parameter in the incoming state is the phase 0o — <f>m- The outgoing state may be completely 
characterized by a similar phase <f) ut, and by the energy ratio rj = E in /E out , both of which 
can be expressed as a function of the incoming phase. Figure E] shows all the information 
needed to describe the classical transition, i.e., the outgoing energy and phase in term of 
similar ingoing quantities. In the case of the energy ratio r\, one can determine the energy 
production for an incoming mode with fixed amplitude but unknown phase, or equivalently, 
the average energy production weighted by the Liouville measure. By integrating the curve 
in Figure El we find 

<V>=^~ f\dfr n = 2.12. (36) 

Z7T JO 

Before moving into the regime close to the singularity let us briefly comment on the 
consequences of including a non-zero center of mass velocity. As we increase the initial 
center of mass velocity v z , the energy production decreases so that the plot for r] (Figure EJ) 
shrinks in the vertical direction. This is to be expected since the z component of the string's 
velocity means that it is closer to the speed of light and hence suffers less of a 'kick' as it 
passes through t — 0. This translates into less energy production or energy loss. 



15 



-jty_ 

-K 



I -Jt 



FIG. 5: The outgoing phase (f> ou t measured at large times is shown as a function of the incoming 
phase 4>i n (left plot). There are two different regions, one where 4> ou t varies steeply and the other 
where it is almost flat. The right hand plot shows the energy production ratio rj = E out /Ei n as 
a function of the incoming phase The ratio r\ is less than one in the steep region of the left 
hand plot, and greater than one outside it. Two neighboring regions are related by time-reversal 
(which sends 77 to 1/??), and they are separated by special points at which 77 = 1. These special 
points, denoted by the vertical dashed lines, correspond to the loop crossing t = with either zero 
momentum or zero radius. 

B. Behavior near the singularity 

Classically, the state of a circular loop with zero center of mass velocity is specified by 
two numbers, its radius and the canonically conjugate momentum. However because of the 
scaling symmetry discussed above, there is really only one physical parameter. Near the 
singularity, we can choose the scale-invariant combination \po/pl\, where p and p are the 
values of p(t) and p(t) at t — 0. This combination compares the radial momentum to the 
size of the loop at the singularity: if \po/pl\ <^ 1) then the radius is large and not much 
changed during the transition across the singularity. In contrast, if |po/Pol ^ 1> then the 
momentum is hardly changed during the transition. In each case, the conjugate variable 
undergoes a large change across the transition. Heuristically, one can think of the transition 
as 'measuring' the loop radius or its momentum in the two cases, so that the conjugate 
variable acquires a large jump, as a consequence of Liouville's theorem. Figures 6 and 7 
illustrate the two situations. The former situation, where the radius is 'frozen' during the 
transition, describes the rising portion of the r)(<f>i n ) curve in Figure 5. As can be seen from 
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FIG. 6: Both numerical simulations and analytical approximations for the radius p (left plot) and 
the canonical momentum (right plot) are shown as a function of time for a circular loop. Due to the 
small invariant parameter |j?o/Pol = 0-l> ^ is the loop radius (left plot), not the momentum (right 
plot) that 'freezes' during the transition. The expansion in the string tension (see text) works well 
in this case, even beyond Hubble crossing. In contrast, the free particle approximation only agrees 
with the numerics close to t = 0, where there is drastic change in the string momentum. If the 
loop radius is taken to be of order the string length l s at a time t s , then this case corresponds to 
6» ~ 1/260. 

the plot, this is the more common situation for initial states with uniformly chosen <fi in . 

Likewise, certain features of Figure 5 are readily understood. Recalling that the phase <f) 
and the 'energy' V denned in (JH2J) are canonically conjugate, one sees how the squeezing of 
the outgoing phase causes the loop energy to be amplified, and vice versa. 

C. Analytical treatment: an expansion in the string tension 

To understand the behavior near the singularity better, we consider expanding about 
t — 0. The simplest possibility is merely to perform a Taylor series in t, but this expansion 
is poorly convergent. Instead, we have developed an analytic solution method which is 
formally an expansion in the string tension, or 1/a'. 

To perform this expansion, let us introduce a formal parameter A multiplying the tension 
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FIG. 7: Another circular loop example, where in contrast to the situation in Figure 6, the ratio 
|po/Pol = 10 so that the momentum p is nearly frozen across the transition, whereas the radius p 
changes significantly - undergoing a complete oscillation. In this case, the expansion in the string 
tension works only until Hubble crossing, and describes the maximum and minimum with less 
precision. If the loop radius is taken to be of order the string length l s at a time t S: then this case 
corresponds to 6q ~ 1/128. 



term in the Hamiltonian, so that 



H — > yp 2 + Xt 2 p 2 . 



(37) 



We now calculate the Hamiltonian equations as before, using (|37|) to eliminate the A 
dependence of the equation for the time-dependence of the Hamiltonian. We obtain 

dH 2 



dp t 2 p dp p 

di = ~ IT' ~dt = ll ) 



dt 



- t (H 2 -p 2 ). 



(38) 



In this form, A appears in only one term in the equations of motion. We can solve the 
equations as a power series in A by inserting the lower order solution into the term involving 
A at each new order. At the end of the calculation, we set A = 1 and fix all remaining 
integration constants. Notice that we do not allow powers of A in the initial conditions: the 
approximation applies only to the dynamical evolution. 

The zeroth order solution is found by setting A = in equations The first equation 
implies p = po, a constant. The last equation then implies H = \Jpl + C 2 t 2 , with C a 
constant. Finally, the equation for p is easily integrated. Finally, we set A = 1 and identify 
the integration constant C . Comparing the expression for H in our solution with the exact 
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Hamiltonian, we can see that C = po, the radius of the loop at t = 0. Hence we obtain the 
zeroth order solution 




This solution is the same as that for a winding string, behaving as a particle of mass M oc \t\, 
which reaches the speed of light instantaneously at t = as described in Ref. 0. 

It is straightforward to compute the solution to higher order in A, although the integrals 
become increasingly difficult. Details are given in Appendix El In general, the series solution 
for p and p to a given order are truncated polynomials in po/Po, where the coefficients are 
functions of x = sinh -1 (u^)- The series converges for much larger times than a simpler 
Taylor series in t in part because x only grows logarithmically with time. If p consists of 
a polynomial of order n in po/Po, then p consists of terms up to order n + 1. Therefore, 
if \po/Po\ 1 the higher power terms will be more important than the lower power ones, 
and vice versa for \po/ p%\ <C 1. To second or higher orders in the expansion, the time- 
dependent series coefficients include polylogarithm functions, which appear naturally in 
Feynman diagrams jsjj, and in number theory. This may be a sign of deeper underlying 
simplicity. 

In Figures 6 and 7, the expansion in the string tension, taken to second order, is com- 
pared with numerical solutions. The analytic approximation shows good agreement with 
the numerics up to times where the loop is starting to be well described by the appropriate 
adiabatic flat-space oscillatory solutions, i.e., when the loop is inside the Hubble radius. We 
conclude that the expansion in the string tension which we have defined is a powerful tool 
for studying classical evolution right across the transition. 



IV. ROTOR 

The classical analogue of the graviton is a rotor solution, a spinning doubled line whose 
ends move at the speed of light. The flat space-time solution representing this incoming 
state is given by equation or, if it has nonzero center of mass velocity, a Lorentz- 
boosted version of it. The appropriate initial conditions for the time-dependent background 
we study are given by computing the flat space-time left and right movers, and translating 
these into expanding universe left and right movers at the initial time. The solutions take 
the form of a doubled line for all time and the two end-points always move at the speed of 
light. The evolution is more complicated than that for a circular loop because there is a 
non-trivial dependence on a. 

There is a set of measure zero on the space of initial conditions where the evolution is 
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FIG. 8: An incoming flat space-time rotor solution (|23|) . with a small initial center of mass velocity, 
v z = 0.01 and t s = 25. The rotor evolves into an "S" shape in the XY plane as the two opposite 
arms speed up towards the speed of light. The outgoing state involves a superposition of higher 
oscillation modes of the string, as may be seen from its evolution. 

singular, when the center of mass momentum of the rotor is precisely zero. Since there is no 
preferred direction orthogonal to the plane of the rotor, its central point must remain at rest 
for all t. This conflicts with the requirement that all points must move at the speed of light 
at t — 0, hence the solution must go singular. However, we do not believe this will cause 
any problem in the quantum theory. States of zero momentum are of zero measure on phase 
space. So we shall study the behavior of a rotor with nonzero center of mass momentum 
p z , in the limit as p z tends to zero. We shall find that the resulting non-analyticity in 
the classical solution is rather mild, and thus likely to overwhelmed by the measure in any 
physically realistic calculation. 

By boosting the static solution (|23[) we can obtain a solution with arbitrary v z , as defined 
in ([25)1 . Any simulation is then characterized by the starting time t s and the initial center of 
mass velocity v z . Contrary to the circular loop, the oscillation phase of the incoming rotor 
is of no physical significance since it can be removed by a spatial rotation. Therefore, as far 
as the classical dynamics is concerned, the solution for a rotor depends only upon v z and 
there are no other parameters to consider. 

The evolution of such a moving rotor across t — is illustrated in Figures |H1 and M In 
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FIG. 9: The profile of the incoming rotor shown in Figure |SJ in the z direction. The passage 
through t = results in a kink being produced at t = 0, which then propagates outwards across 
the string. 



the XY plane, it develops an ! S' shape. This may be understood as a consequence of the 
opposite arms of the rotor speeding up to the speed of light in opposite directions as t = 
approaches. In the Z-direction, as the central point speeds up to the speed of light it creates 
a kink which then runs out across the string. 



A. Profile of the Kink 

As the rotor approaches t = 0, it develops a kink in the direction of its motion, whose 
size and shape depends on the magnitude of the transverse momentum. For zero transverse 
momentum, the classical evolution is ill-defined. This is a set of measure zero, but the break- 
down of the classical equations at p z = may be indicative of some divergence. Therefore 
it is important to study the evolution for small p z to see whether physical quantities diverge 
that limit. 

As the transverse momentum gets smaller, the kink gets narrower and narrower, as shown 
in Figure ITU1 

The profile of the kink may be modeled analytically as follows. Consider the equation 
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FIG. 10: The kink in the transverse velocity z around the middle of the rotor, for two very small 
velocities. The profile Q41JI fits the shape very well for small transverse velocities (red dashed line), 
where the quantity < Zq >= J e^z^do j J e^da measures the center of mass velocity at t = in the 
z direction. The initial conditions were: t s = 5, v z = 0.05 for the wider curve and v z = 0.01 for 
the narrower curve. 



for the velocity of the string in the z direction, z. From equation ([10)1 . at t = we have 

vr,(0) ir z (Q) 



z = 



6(0) 



7T 



(0) : 



(40) 



where n = ex is the canonical momentum density at each point on the string, and we used 
fTTj) to get the last equality. For small transverse velocities, near the center of the rotor we 
can replace tt z ~< z >= / e^z^da / J e da and use the Taylor expansion 7r^ + 7r^ ~ (C a) 2 , 
with C a constant. Therefore, near the center of the rotor, the kink's profile at t = is 
approximated by 

z = . 1 (41) 

This fit works well against our numerical results, as shown in Figure with C 2 ~ 0.2. 
As we take the centre of mass momentum of the loop to zero, the kink becomes more 
and more strongly localized and involves higher and higher oscillation modes of the string. 
However, these do not contribute significantly to the energy. Figure [TT] shows r\ = E ou t/E in 
as a function of the loop momentum. In the limit of small momentum, rj tends to a finite 
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FIG. 11: The energy ratio r/ = E ou t/Ei n for a rotor is plotted against the average center of mass 
speed of the loop at t = 0, i.e., < zq >= J e^z^da / J e^da, showing that the net energy production 
is unaffected by the non-analytic spike produced in the limit of small transverse momentum. 

constant ~ 1.93, showing that there is little energy associated with the spike generated at 
the center of the rotor. We conclude that at low momenta, we produce a a non-differentiable 
but finite spike, although physical quantities like the energy or momentum remain perfectly 
finite. 

B. Expansion in the String Tension 

Just as we have done for the circular loop, we can describe more general string states, 
including the rotor, using a formal expansion in the string tension. The details are given in 
Appendix [XJ Following the same general method as given for the circular loop, the zeroth 
order solution is given as 



where p , x and x' Q are functions of a. 

Higher terms in the expansion are more difficult to calculate and the integrals can only 
be done numerically (see Appendix El for details). Figure [T21 shows the numerical solution 
plotted against the zeroth order term in the expansion, the first order term, and the flat 
space approximation, respectively. The comparison is made for an end point of the rotor 
(right plot) and a representative point further in (left plot). For the former, the adiabatic 




(42) 
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FIG. 12: The numerical solution for two points on a rotor. The left plot shows the x coordinate at 
a = 7r/4 and the right plot shows the x coordinate of the end point a = 0. The adiabatic flat space 
approximation reproduces the numerical solution reasonably well until the loop crosses the Hubble 
radius. Beyond this point, an expansion in the string tension becomes a much better description. 
In this plot, the initial center of mass speed was v z = 0.05 and t s was taken as 35. 



approximation holds very accurately up to Hubble crossing. Since the point is moving at 
the speed of light, the zeroth and first order terms of the expansion are nearly identical, 
providing a good approximation to the motion outside the Hubble radius. At Hubble re- 
entry, the solution reverts to a form close to flat-space evolution but this is more complex 
to determine hence we have not attempted to graph it. For generic points on the rotor, 
the adiabatic approximation is less accurate at earlier times: the speeding up of the central 
point as the kink is created in effect extracts energy from the normal spinning motion. As 
can be seen from the diagram, the first order expansion in the string tension does reasonably 
well in modeling the behavior around t = up to the first turning point. 



C. Outgoing state 

When the rotor is evolved forward to large times, the outgoing state consists of a com- 
plicated mixture of higher oscillation modes. The precise mixture is most easily understood 
by representing the outgoing solution in left and right moving modes, which tend to fixed 
curves on the unit sphere at late times. One can also use the right and left movers to under- 
stand the motion of the rotor across the singularity (Figure The two curves start out 
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FIG. 13: The evolution of the left mover I for an incoming rotor. The left mover defines a closed 
curve on a three-sphere. The right mover follows the same curve but with the opposite orientation. 
The condition Q15|) that the left and right movers coincide t = 0, means that the evolution must 
deform the trajectories over the north pole of the three sphere. The only way to get the two curves 
to coincide is by deforming them into a double line at t = 0, running directly through the north 
pole (or the south pole if v z is negative). The initial parameters here were v z = 0.05, t s = 45. 

near the equator (for small v z ) and oppositely oriented. For zero center of mass velocity, the 
the curves are confined to the xy plane and there is no way for them to coincide at t = 0. 
However, if v z is positive (or negative), the curves sweep over the north (or south) pole and 
coincide at t — 0. 

To quantify the level of excitation produced in the passage across t = 0, we can track the 
the evolution of the right and left movers in time at one particular value of a on the string. 
For very large times, the right and left movers look like flat space solutions and become 
periodic functions of proper time. Once the loop is well inside the Hubble radius, for a 
single period the difference between conformal time and proper time is negligible. Hence we 
can just choose one value of a and follow the evolution of left and right movers there. We 
write the expansion of the left and right movers (|21j) as 

+oo +oo 

r *= a^e-^-^ , f = ]T o&e"* 1 ^, (43) 

n=— oo n=— oo 

and compute the a l n and a l n by Fourier transforming the solution with respect to t over a 
single period. As is seen in Figure 19, identifying the periodicity is straightforward. As the 
initial center of mass velocity v z is tuned down, modes at higher and higher n are excited, 
although as shown above, the total energy in the outgoing string converges to a finite limit 
as v z tends to zero. 
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FIG. 14: Fourier amplitudes of outgoing right movers calculated in the asymptotic adiabatic regime 
for different initial transverse speeds. The red dashed line shows the z component only, hence a$ 
(the n = value of the red dashed line) is the final speed v z . By symmetry, a x and a y have only 
odd-n modes whereas a z has both odd and even n components. In these plots t s = 25. 

V. CLASSICAL ANALOG OF THE ANTISYMMETRIC TENSOR 

To complete our discussion of low mode solutions, we consider configurations with in- 
termediate angular momentum, which are the analog of the massless antisymmetric tensor 
states of the quantized string. The corresponding flat-space solution was given in (|24)l. and 
its evolution, illustrated in Figure 2, consists of a spinning state whose shape alternates 
between an ellipse and a doubled line. The motion is simplest when projected onto planes 
spanned by the vectors: x = (1, 0, 0), v = -^(0, 1, 1, 0, 0) and w = ^(0, —1, 1, 0, 0). 
In the xv plane the solution appears as an ellipse which shrinks to a point and re-expands. In 
the xw plane the motion is similar to that of the rotor, but the length of the rotor oscillates, 
(see Figure EI). 

If we follow the evolution of such a solution towards t = 0, once it has crossed the 
Hubble horizon, it practically freezes in size in the x and w axes, and either its position or 
its momentum freeze in the v coordinate. Apart from the special case where the elliptical 
motion has a maximum at t — (which is mildly non-analytic, like the case of a circular 
loop), the motion is regular across t = and produces a tower of higher excited modes (see 
Figure ITUj) . If the position is the frozen variable in the v direction, there are fewer excited 
modes produced, and conversely if the momentum is the frozen variable. As for the rotor, 
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FIG. 15: The energy ratio r\ for each mode number n, representing the amplitude of each outgoing 
mode divided by the amplitude of the incoming mode. Larger initial speeds v z generate smaller 
amplitudes in the higher modes. The decay of rj with n resembles the plane wave estimate (~ 
e^v 7 ™ i 'y/n) found in [5j for high velocities, but at low velocities rj decays more rapidly than 
expected from that work, and in particular the production of high n modes remains finite. 




X JC j£ jC jC jC 

FIG. 16: Evolution of the classical analog of the antisymmetric tensor 1)24(1 across the singularity. 
The six plots on the left show the behavior in the xv plane, where, inside the Hubble radius, the 
solution is an oscillating elliptical loop. Higher modes are excited as t = is crossed. The plots on 
the left show the projection into the xw plane, in which the solution always appears as a doubled 
line. Again, it is clear from the pictures that higher modes have been excited. For this case, t s = 6 
was used. 
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FIG. 17: An incoming classical analog of the antisymmetric tensor state is represented by left 
movers (blue lines) and right movers (black lines) on the sphere. At t = the left and right movers 
coincide in a doubled line. 

we can picture the evolution of the antisymmetric tensor's classical analog in terms of left 
and right movers on the sphere (Figure IT7j). 
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VI. QUANTUM EVOLUTION ACROSS T=0 



In this section we treat the quantum production of excited string modes in a simple 
approximation, which may be extended to discuss more complicated string states. As we 
have already indicated, for large string states it is reasonable to expect that, provided the 
quantum theory remains well-defined, quantum corrections should be small and our classical 
calculations should be accurate. 

Here we shall limit ourselves to an illustrative example, consisting of a truncation of the 
full theory to the simplest oscillation mode, namely a circular loop. The method we use 
generalizes to more interesting excited states and we shall study this in future work. We 
start from the classical equations of motion for a loop, 

dp f_p dp p_ 

dt ' dt H' P + P ' 1 J 

We would like to quantize these equations, but the nonlinearity renders this difficult. We 
proceed in the following approximation. The operator H 2 = p 2 +t 2 p 2 is positive definite, and 
so it is not unreasonable to attempt to approximate the terms H in the first two equations 
with the c-number quantity, 



H c (t) = ^J(H 2 (t)). (45) 

With this replacement, we have only to solve the linear quantum system whose Heisenberg 
equations of motion are 

dp = t 2 p dp = p 

dt H c (ty dt H c (ty 1 ' 

where H c (t) is as yet undetermined. 

The general solution to (J44j) may be expressed in terms of two linearly independent 
solutions, defined at some initial time to, as follows: 

p(t) = Pofp(t,to)+Po9 P (t,t ), 

P(t) = Pof P (t,t ) + Pog p (t,t ), (47) 

where f p , f p , g p and g p satisfy the initial conditions f p — 1, f p — 0, g p — and g p — 1 at 
time t = t . 

We now wish to express the operators po an d po in terms of creation and annihilation 
operators for the incoming vacuum state. To do so, we need to identify the relation between 
p and p and the properly normalized coordinate and momentum of a harmonic oscillator at 
large \t\. Noticing that the quantity as 

2t 2t ' v ; 
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FIG. 18: The quantum creation of higher excitations on a circular loop, sent in in its ground state. 
The upper plot shows the two independent solutions f p and g p defined in the text. The lower plot 
shows the evolution of the 'energy' H c /(2\t\). This is not well-defined at t = 0, even though the 
equations of motion are regular there. For a loop sent in in its vacuum state, with an energy | in 
these units, the outgoing state has an energy of approximately 2.54. The increase in energy is due 
to quantum mode creation as the loop crosses t = 0. 



defined in (J32|) . is an action variable at large \t\ (see Section UlI A|) . we identify the proper 
momentum 7r = p/ y\t\ and radius R = p\J\t\ as the asymptotic coordinate and momentum 
of a corresponding harmonic oscillator. These variables can then be expressed in terms of 
creation and annihilation operators, denned at some large negative time t : 

(a + %) 



Po 



Po 



\to\ (do 



(49) 



'2|tb| V 2 

It is then clear how to compute the function H c (t), from the formulae above, for any chosen 
incoming state. In particular, in the incoming vacuum state defined by ao|0)m = 0, we find 

1 / fn , . , O 1 ( j ~g 



m) = 2[h +ltolg2p+t \\to 



9 + 1*0^ 



(50) 



It is straightforward to solve equations (|4*Cjj) for the two independent solutions / and g, using 
(150 j) to define H c (t). This then provides a self-consistent semi-classical approximation to the 
full quantum evolution. 

Figure shows the results of this calculation, for the incoming vacuum state. From 
formula ()48|). we identify H%/(2\t\) as the asymptotic harmonic oscillator energy at late 
times, which may in our chosen units be expressed as i n (0\(N out + |)|0)j„ where iV is the 
usual number operator. Therefore, from the Figure, we see that in the outgoing state, 



1 (0|iVo,t|0) m = |/3| 2 ^2.04, 



(51) 
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where (3 is the usual Bogoliubov coefficient measuring the amount of particle creation [16]. 
Similar results are obtained for excited incoming states. One can straightforwardly compute 
the Bogoliubov coefficients a and (3 including their phases, using the same method. The 
results will be given elsewhere. 

VII. SUMMARY 

We have herein discussed how classical strings propagate through the simplest possible 
big crunch/big bang transition in M theory, and made a start on considering the quantum 
problem. As we have discussed, the usual a' expansion must be replaced by an expansion in 
the string tension, i.e. I /a', once the Hubble radius falls below the string length. We have 
presented an analytic solution of the Nambu equations as a formal expansion in the string 
tension, which could form the basis for a quantum treatment of the small time regime. 

Using the decomposition into right and left movers, we have shown how higher oscillation 
modes are excited by a finite amount as the string crosses singularity. We have also given a 
discussion of some special cases where the string is constrained by symmetry to be static at 
t = 0. Except for a set of configurations of measure zero, the classical counterparts of the 
dilaton, graviton and antisymmetric tensor string states evolve smoothly across t = 0. 

We have made a modest start at describing the quantum evolution of string across t — 0, 
discussing a truncation of the theory to the lowest mode describing a circular loop. The 
method we have employed appears to be readily extendible to more complex incoming states, 
and to the inclusion of more and more string modes. As this is done, we can begin to address 
the subtle issues of renormalization which we have so far ignored. It may also be feasible to 
go beyond the first semiclassical calculation we have described. 

Our main conclusion is that classical strings pass smoothly through singularities of the 
simple type studied here, and the excitations they acquire in doing so may be readily calcu- 
lated in a formal expansion in the string tension. While much remains to be done to develop 
a full quantum description, our findings reported here are encouraging. 
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APPENDIX A: SERIES EXPANSION AROUND t = 



1. Circular loop 

As explained in Section IIII CI in order to study the motion of the string across t = 
we perform a formal expansion in the string tension. We do so by introducing a formal 
parameter A in front of the relevant term in the Hamiltonian, and then we solve the equations 
of motion to a given order in A. Finally, we take A = 1. With this parameter A, equations 
(J2HJ)-(|IID) take the form (jHHj) . which we repeat here: 

dp _ t 2 p dp _ p dH 2 _ 2 2 2 

dt ~ - x ir> &-h> ~ i {H ~ v ] - (A1) 

From this point, there is only one consistent way to proceed. At lowest order in A, we solve 
the above equations with A = 0. Then we substitute this solution into the term involving 
A and integrate with respect to t to obtain first p, then H and finally p to order A. This 
solution is again substituted into the right hand side of the first equation, and the procedure 
continued. Effectively, we have 

dpM _ At2 /M(-> d(r 2 (H 2 )^) _ 2 ( P 2 )^ (A2) 



dt \HJ ' dt t 3 

which can be written as integral equations, 

p(») = Po + A/ (n - 1} , (H 2 )^ =p 2 + p 2 t 2 - AXpoJ^-V - 2X 2 K {n - 2 \ (A3) 

where p(0) = p , p(0) = p and I {n \t), J( (n) t), and K^(t) are defined as 



(A4) 



wo.f («JU (A5) 

Jo t 6 

Now equations ()A3j) have the desired form. Finally, we use these solutions for p and H to 
construct p to the same order in A by expanding the quotient 

P I P A . (A6) 



dt \H 
To zeroth order in A, we find 



p<®=Po, H^ = ^p 2 + CH 2 , (A7) 

where C is an integration constant that we will fix at the end, after setting A = 1. The 
solution to p(°) is then given by 

^^^^'Sl (A8) 
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In order to construct the next order solutions, there is a change of variables which simplifies 
the expressions, given by 

X = smh- 1 ( f\t) , (A9) 



bo I 



thus the zeroth order solutions look like 



p<°> = po, # (0) = Po cosh( X ), = po + ^X- (A10) 



To next order in A, we get 



P (1) = Po 



(# 2 ) (1) = vl 



1 + (2 X - sinh(2 X )) + A^j ( X 2 + sinh 2 ( X ) - X sinh(2 X ) 

cosh 2 (x) + A^^(2xcosh(2x) - sinh(2x)) 



2|C| 3 

+A^(x 2 cosh(2x) + sinh 2 (x) - xsinh(2x) 



(All) 



The parameter A governs the dynamics of solutions near t — 0, so it can not appear in the 
initial conditions. Consequently, after setting A = 1 the value for C is determined using the 
equation for the Hamiltonian (|29p. which implies H"(0) = 2p\. From the solution of (H 2 )(°\ 
or equivalently (H 2 )^\ at t — one reads C = po- Therefore, the solutions for p^ and 
are just polynomials in the scale-invariant parameter |p / Po I • 

Note that the functions p^ has an odd and an even piece with respect to x- Therefore, 
the function will be more symmetric around t = if |po/Pol ^ 1> an< ^ more asymmetric for 
the converse value. The odd-x part will always have a zero near t = 0, which represents 
either the maximum or the minimum in p near t = 0. The other zero of p depends on 
the even-x piece, which is close to t — only for |po/Pol ^ !• However, the series is only 
accurate for small \ which guarantees —t 2 p/H is small after taking A = 1, so only if the 
invariant parameter |po/Pol ^ 1 one can §° a little beyond the Hubble radius crossing. 

By expanding the quotient on the RHS of equation (jA6|) to first order in A we obtain 
which looks like 

i 

po\cr ' "4c 4 



P (1) = Po 



1 + -^rX + ( x(2 tanh(x) - x) - smh 2 (x)) (A12) 



3 

A - P ° (9 X + 2 X 3 + 3 X cosh(2 X ) - 3 sinh(2 X ) - 6 tanh(x) - 6 X 2 tanh(x) 



12p Q |C| s 

Again, has piece which is an even function of x-, Du t now there are two terms which 
are odd-x functions. It is due to this extra order term in po/p 2 , that the p solution is more 
symmetric if |po/Pol ^ 1 an d asymmetric if \po/Po\ ^> 1, in opposition to p^\ 

In order to obtain better precision and more maxima and minima in p with this series, 
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one has to go beyond first order in A. However, the integrals become harder and harder to 
solve. Just to illustrate the flavor of the next order solutions, is given by 



P 



(2) 



p cosh 3 (x)< cosh 3 (x) + -4 [xcosh 3 (x) 
I Po 



(A13) 



+ 



Po 



cosh 2 (x) 



+ Ax 2 ) cosh(x) + cosh(3x) - 8xsinh(x) 



_ ( .^ (3 [32cosh 3 ( X ) (tt 2 + 12 X 2 - 24 X ln(2cosh( X ) + 12Li 2 (-e~ 2 

+X(9 cosh(x) - 4 cosh(3x) - 12 cosh(5x)) + (24^ 2 - 42) sinh(x) 
+ (120x 2 - 27) sinh(3x) + 15sinh(5x) 



+ 



Po 



3 (347 + 208x 2 ) cosh(5x) + 3 cosh(7x) - 1440x sinh(5x) 



12288p 7 L 

+ 128 cosh 3 (x) (l44x 2 ln(2) - 96x 3 + 5^ 4 + 48(1 + 3^ 2 ) ln(cosh(x)) 
+6 cosh(x) (-647 - 40^ 2 + (473 + 1248) cosh(2x) 
-96 X (2(17 + 6 X 2 ) sinh(x) + (65 + 28 X 2 ) sinh(3 X ) 
-96 (24cosh(x) 3 (4Li 3 (-e~ 2x ) + 8Li 2 (- e - 2 *) + 3((3) 

4cosh(x) (l6cosh 2 (x)(40vr 2 + 7tt 4 + 480x 2 + 720^ 4 - 20x 

-320x 3 cosh(x) (-4 + 13cosh(4x) + 57cosh(2x) + 192cosh 2 (x) ln(2 cosh(x))) 
+ 15360 cosh 3 ( X ) ((2 + 6x 2 )Li 2 (-e- 2x ) + 6xLi 3 (-e" 2x ) + 3Li 4 (- e - 2 *)) 
+30 sinh(x) (l204 + 1824^ 2 + 320x 4 + 1845 cosh(2x) + 2560x 2 cosh(2x)) 

+30 sinh(x) (448% 4 cosh(2x) + 516 cosh(4x) + 480x 2 cosh(4x) + 3 cosh(6x) 



Po 



122880pg 



where we have already set A = 1 and C = po- in the above expression, ((x) is the Riemann 
Zeta function and Li n (z) is the n-degree Polylogarithm function, which is defined as [2^ 



U n( Z ) = 



00 z k 



(A14) 



k=l 



where z belongs to the open unit disk in the complex plane, and is defined uniquely for 
\z\ > 1 by analytic continuation. 



2. Expansion in the string tension for a general string configuration 

In this Appendix we want to construct a formal expansion in the string tension for an 
arbitrary string configuration. Our starting point is the Hamiltonian corresponding to the 
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action (JUJ), 

H = J doH, H = ^V 2 + Xt 2 {x') 2 , (A15) 

where in analogy with the circle, we have chosen to work in units where H2O0 — 1 but we 
have inserted a formal parameter A in front of the term involving the string tension. The 
Hamiltonian equations read as follows: 

d t v = \t 2 d a (^pj, d t x = ^, d t n 2 = - t {n 2 -v 2 ). (Ai6) 

The second equation can be written in the integral form 

x = n , (A17) 
where po = Po{cf) is the momentum distribution at t — 0, and 

^•^/.Mwf (A18) 

Similarly, TC(t, a) also has an expansion with respect to A, given by 

(H 2 ) {n) = po 2 + D 2 t 2 - AXp ■ J^- 1} - 2XK {n - 2 \ (A19) 
where D 2 , is an integration constant that we will fix after taking A = 1, and 

J {n) M = £p n) <?>°)> K{n) M = jf |([^ ff )] 2 ) w ( A2 °) 

With these definitions, equation (|A17|) becomes d t x^ = (j^J \ and hence the zeroth order 
equation in A is 

iM = _ 



S m = , P ° , (A21) 



which is easily solved by 

f (o) = f + A- sinh- 1 ( ^t] , (A22) 
\D \ \ \po\ J 

where xq is the string shape at t — 0. After setting A = 1 the Hamiltonian equation IA15I 
implies H e "(0) = 2(x' (0)) 2 , which translates into D = \x'(0)\. 

Now, to first order in A, we find the equation for x is given by 



1,(1) Po 
x 



Po + Dlt 2 



1 + A/<°> + A 



2t 2 p ■ J(°) 



Po 2 + D 2 Q t 2 



(A23) 



and with a similar change of variables to that we introduced before for the circular loop, 
X = sinh" 1 (^p)> the first order solution reduces to 
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These integrals can be reduced to a single integral by changing the order of the integrals, 
provided no convergence issues arise. For situations where the parameter Ipol/^'ol ~ T there 
is no difficulty and the the integrals can be exchanged in the following way 



[ X ''dtil® (•&,&) = -4 / X ^sinh 2 (^) cosh(^)(^ - $)d a ( ?ZE^M 

Jo -Do Jo \ H($, a) 

coshW^ ^ - _ f s . r sinh 2 ^)] a (d a x& a) 

Jo sinh 3 ^ 1 ' ] ~ 2Dlk mC0Sh ^ 



sinh 2 ^)' 
sinh 2 (-$) 



(0) 



(o) 



We have used these last expressions for the example in Figure (jl2j) . 

In principle, the method we have presented may be extended to second or higher order 
in the string tension. However, the integrals become very complicated and only in simple 
cases, like the circular loop, one could actually calculate them analytically. 



APPENDIX B: CHECK OF NUMERICAL EVOLUTION 



In order to check our numerical algorithm for following the string evolution, we have 
performed a number of tests. Figure shows the behavior of right and left movers for 
an incoming rotor. One can readily distinguish the regime of adiabatic oscillations from 
the super-Hubble evolution around the singularity. In the former region, the right and left 
movers oscillate as in flat space-time. Different positions o on the string are distinguished 
only by a phase difference. Moreover, because the rotor is isotropic initially, the incoming x 
and y components of the right and left movers are equivalent. However, after the singularity 
these components behave differently, in agreement with the loss of isotropy in the xy plane. 
Furthermore, the spinning doubled line structure makes the left movers equivalent to the 
right movers, with only a phase difference. This structure is preserved across the singularity, 
hence only a n , with n-odd, are produced during the transition. Contrary, the z direction 
gets modified during the singularity, from constant right and left movers to a mix of excited 
states with both, odd and even n. In the 1/a 1 region, the evolution of right and left movers 
is completely non-linear and there is not any correlation between right and left movers, 
between different components and, perhaps more interesting, between different points in the 
string. 
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